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Abstract. We explicitly calculate the moments i„ of general Heisenberg Hamiltoni- 
ans up to sixth order. They have the form of finite sums of products of two factors, the 
first factor being represented by a multigraph and the second factor being a polynomial 
in the variable s(s + 1), where s denotes the individual spin quantum number. As an 
application we determine the corresponding coefficients of the expansion of the free 
energy and the zero field susceptibility in powers of the inverse temperature. These 
coefficients can be written in a form which makes explicit their extensive character. 



1. Introduction 

The moments of the general Heisenberg Hamiltonian describing a large class of spin 
systems represent valuable information for various purposes. They can be used to check 
numerical exact or approximate calculations of the energy spectrum of such systems. 
Further they directly appear in the high temperature expansion (HTE) of the partition 
function. Hence they can be used to calculate the leading terms of the HTE of physical 
relevant functions such as specific heat or magnetic susceptibility, see e. g. [1] or [2]. 
Analytic expression for the moments in the general Heisenberg case (arbitrary number 
of spins and coupling constants J^^, arbitrary spin quantum number s) are interesting 
since they can be used to investigate the dependence of physical properties on the 
coupling constants even in cases where a numerical investigation is not possible due 
to large N or s. It seems that such analytical expressions for moments have only be 
published up to order three, see [3], although the polynomials which appear in such 
expressions are known up to 8th order, see [1]. Those papers which consider higher 
order expansions are usually confined to special coupling geometries or and/or s = 1/2, 
see e. g. [T7]. 

In this article we only give the analytical results for the moments t^, t^, without 
explaining the method by which we obtained them. This will be done in another pa- 
per. From these results we have determined the coefficients of the susceptibility's HTE 
C4) C5, Cq by a method which is briefly sketched. Similarly, the coefficients a„ of the free 
energy's series expansion (or, more precisely, of —f3F{f3)) are calculated up to 6th order 
in /3. We have checked our results by comparison with numerical results for spin systems 
with random coupling coefficients and N < 7 and s = 1/2, 1. Moreover, our expression 
for C4 has been checked by comparison with the published results for chains [18] and 
Ji — J2 square lattices [T9] [20]. 
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We consider the following notations: 



N 

Sn 

s = 
r = s{s + 1) 



1 1 3 

2 ' ' 2 ' 



Jfiu = Jiyfj., I < fJ' 7^ < N 



E 



1,2,3 

N 



H = E 
(2s + 1) 

, _ Tr(H") 
~ (2i+Tj^ 
_ Tr(S(3)2H") 
Z^" — (2s+l)^ 

13 = ^ 

UR\ - o Tr(S(^)^exp(-/3g)) 
A\P) P Tr(cxp(-/3H)) 

XW) = EZl CnP- 
-(3F{(3) = In (Tr e"^^) 

-l3m = En=o^nf3- 



Number of spins 

Group of permutations vr : {1,2,..., A^} — )• {1, 2, 

Single spin quantum number 

Abbreviation 

Coupling constants 

Spin vector operator of the /ith spin 

Total spin vector 

ith component of the total spin vector 
Heisenberg Hamilton operator 
Dimension of the Hilbert space 
Normalized moments of H 

Normalized magnetic moments of H 

Inverse temperature 

Zero field susceptibility 

High temperature expansion of xiP) 

Free energy F{P) 

High temperature expansion of F{(]) 



In order to write our results in a compact way we will utilize some graph theoretic 
notations. This is a common practice when dealing with series expansions, see e. g. [2]. 
Let ^ be a multigraph consisting of g nodes (vertices) and J\f{i,j) = Af{j,i) bonds 
(edges) between the ith and the jth node. We do not consider "loops", i. e. M{i,i) = 
for al\i = 1, . . . ,g. The total number of all bonds, 7(^) = Ej<i -^l^' "^^^^ called the 
size of Q. Q is not necessarily connected, see the examples of multigraphs below. We will 
identify the set of g nodes with {1,2,..., g} and the set of spins with {1,2,..., A^}. 
All multigraphs Q^, u = 1, . . . , 86 needed in this paper are represented in tabled 

From the expansion 

= E n J^^^^ Tr ( n • ^'A (1) 

it is clear that the expressions for the moments t„ involve various products of coupling 
constants J^y. The structure of these products can be represented by the multigraphs Q 
defined above, such that the factors J^j, correspond to the bonds of Q with multiplicity 
i. The sum of different products in ([1]) of the same structure will be obtained by an 
evaluation of denoted by for the spin system under consideration. Q denotes a real 
number which depends on the coupling constants and only implicitly on the number A^ 
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of spins. This number will be defined according to the following statements, 
(i) If ^ > we set ^ = 0. 
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(ii) If (yf < we consider the symmetry group G of defined by 

G = {tt G Sn\t^ maps {1, 2, . . . , (?} onto {1, 2, . . . , (7} and 

M{i,j)=N{^{i),'K{j)) for all 1 < I, J < g} (2) 

Further we consider the quotient set S]\f/G consisting of left cosets tt G, tt & and 
choose a representative tt^ from each coset, ne E nG, £ = 1, . . . , L = \Sn/G\ = j^. 

(iii) Then we define 

L 

^ = ^ n i'^^di),Mj)) ''^ • (3) 

e=l l<i<j<g 

Obviously, the definition of Q does not depend on the choice of representatives iri since 
the product ni<i<j<g {'^Tri{i),ne{j))'^^^''''^ is invariant under permutations from the sym- 
metry group vr G G. 

In order to illustrate this definition we consider an example of = 4 spins and 

^ = A , hence g = 3 < A = N. The symmetry group G consists of all permutations 
of {1,2,3,4} leaving 4 fixed, hence |G| = 3! = 6. Thus we have L = || = 4 left cosets 
from which we choose the representatives tti = (1)(2)(3)(4) = id, 112 = (1234), tt-^ = 
(13)(24), 7i4 = (1432). Hence Q = J12J23J13 + J23J34J23 + J34J14J13 + </i4^i2^24- 

The coefficients c„ of the susceptibility's HTE (and similarly the a„ of the free energy 
HTE) will contain products of evaluations Qu Qfj^. These expressions can be simplified 
using rules which transform such products into linear combinations of other evaluations. 

To give an example, we consider ^2 = o = (^YJ^^<u ^i^^^ {Y.k<x ^Ix) ■ 

is obvious that this product can be written as a sum over evaluations of the three 

multigraphs which can be combined from •-• and o ^ namely «=» , *=»— and 

& . Similar expressions can be derived for other products of evaluations where suitable 
coefficients will compensate the change of symmetry groups. For the purpose of the 
present paper we will need the following "product rules": 



Q1Q2 = e?4 + + 6^6, (4) 

Q2 = G9 + 2^12 + 2^17, (5) 

Q2 Qs = Sw + Oi3 + Qu + Gig + ^18, (6) 

Qi Ga =09 + Qio + Qii, (7) 

Q2 Qi = Q23 + ^26 + ^31, (8) 
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GlA + Gl9 + G20, 


(9) 


G2 G4 




G23 + G26 + G31, 


(10) 


Gs Gi 




G24 + G27 + G2S + Gm + ^32, 


(11) 


Ga 




^58 + ^^62 + ^^66, 


(12) 


G2 G5 




^24 + ^26 + 2^33 + 2^34 + ^35 + 2^44 + ^46, 


(13) 


G2 Ge 




^25 + G?,\ + 2^35 + 2^36 + Ga'o + 2^47, 


(14) 


G2 G7 




^28 + ^38 + ^48; 


(15) 


Gs G7 




^34 + ^37 + ^50 + 2^51 + ^52 + ^54, 


(16) 


G2 G7 




^60 + i^64 + 2t/67, 


(17) 


75-2 

G7 




^70 + 2^74 + +2^83 + 2^85, 


(18) 


G2 Gs 




^49 + ^29 + ^30 + ^37 + ^39 + ^42 + ^45, 


(19) 


Gi Gg 




^23 + t/24 + t/25, 


(20) 


G2 G9 




^58 + ^59 + ^61, 


(21) 


G2 G12 




^59 + 3^68 + 3^70 + 2^71 + ^73, 


(22) 


Gi G12 




^26 + ^33 + ^34 + ^35 + ^36, 


/ r\<r\\ 

(23) 


Gi GiA 


— 


^28 + 2^34 + ^37 + ^40 + GaX-, 


(24) 


G2 GlA 




Gm + ^63 + ^69 + ^76 + ^78, 


(25) 


Gi Gi7 




^31 + ^44 + Ga% + ^47, 


(26) 


G2 Gl7 




^61 + G7\ + 2^73 + 3^81, 


(27) 


Gi G21 




^42 + 2^51 + ^53 + ^55, 


(28) 


G2 G21 




^65 + 2^74 + ^77 + ^82- 


(29) 



From these equations one can derive further ones for multiple products, for example: 



G\ G2 = G23 + G2A 

+ 2(^ + ^ + ^ + ^ + ^ + ^ + ^ + ^). (30) 

3. Results 

3.1. Moments 

It turns out that the moments t„ can be written in the following way: 

tn=J2'^- P-^''^ ■ (31) 

Here the Giy, G T„, denote certain multigraphs of size n and the are polynomials 
of order < n in the variable r = s{s + 1). It is crucial that these polynomials depend 
neither on nor on the coupling constants J^^, whereas the terms Gu depend only on the 
coupling constants and only implicitly on N via Q. For the determination of tn it thus 
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suffices to enumerate the multigraphs Q^j, z/ e T„ and the corresponding polynomials 
Py. The first moments are well-known: 

ti = (32) 

^2 = E Ji = ° = ^ V ^^^^ 

The third moment has been calculated in [3] . We reproduce this result using our 
notation and the multigraphs represented in table [H 

^3= -Ir^Gl+lr^g'r. (34) 
3 

The next moments t^, t^, have, to our best knowledge, not yet been published, 
although the polynomials which appear in these expressions are known up to 8th order, 
see [1]. 

They are given by the following expressions: 



r\2 



3r)) g, 



2 



15 



-l + 3r) g 



12 



gi4 



(35) 



and 



K^il + 2(-l + r)r) g23 + ^r%l - 2r) g 

io 



26 



+ V(3 
9 ^ 



20 



28 



9 



r^6?3i + ^r^(-l + 3r)^?38 



27 



10 
— 1 
27 



20 



56 ) 



(36) 



te = r^{32 - 87r + 88r^ - 30r^ + 15r^) g^s + -r^(-3 + 8r - 8r^ + Gr^) g. 

105 6 

- ^r3(l - 3r + 3r2) ^ + lr^(2 - 2r + 3r2) ^ + ^r=^(-2 + 3r) ^ 

+ V(l - 2r) ^ + y r^(l - 2r) ^ + ^r^(3 - 8r + 12r2) ^ 

+ V ^ - ^r^ ^ + ^r^(l - 3r + 3r2) ^ 

+ ^r\l -2r)g^+ lr%-l + 9r - 22r2 + 22r3) ^ 
9 5 

2 2 10 

+ -r\2 - lOr + 15r2) gj, + -/ g^^ + Yr^(-1 + 3r) 6^73 

+ V(5 - 24r + 36r2) ^ - ^ + ^r^(l - 5r) ^ 
+ ^^'(-1 + 3r) ^ - y ^ + ^ - ^ 



Moments of Heisenberg Hamiltonians 10 



^^-r' ^ + ^ + ^r^(-l + 3r) ^3 



9 27 
,24-— 80g-^ 80g-^ 

+ -r ^84 + yr ^§5 + i^r ^86 • (37) 

3.2. Free energy 

It is well-known that the coefficients of the power series for the free energy 

00 

- = In (Tr e"^^) = ^ a„/3" (38) 

n=0 

can be expressed in terms of the moments t„ and its products. As indicated in section 
[21 a variety of product rules can be used to simplify the resulting expressions. This 
simplification, which is sometimes also referred to as the "cumulant expansion", see [2], 
has the further advantage that it reveals the extensive character of the section 
13.31 for a more detailed discussion. The ffist seven coefficients of the series (138|) read as 
follows. 

ao = iV ln(2s + 1), (39) 
ai = 0, (40) 

as = ^, (41) 
a3 = ^r2^-^r3^, (42) 

04 = - —r'^i-l + r + r^) ^ - — ^ - — ^ + — (43) 
4 V ^ ^ 1^9 ^Qg b/14 -r ^21, V ; 



1 „ „ . _ 1 

26 



«5 = - 7t;t7T^^(-3 + 6r + 4r2) ^23 - tttt^^^ ^: 



2160 ' ' 432 

^^r^(-3 + 8r + 8r2) ^ 
1080 ^ ^ * 



+ ^r'fe + 5ijr'a„-ir=fe. (44) 



as = — ^ — r2(192 - 522r - Qlr"^ + 240r3 + 160r^) ^58 



453600 

+ ^^r^(-9 + 12r + ^r"^) + -^r^(-l + 3r + 2r^) ^ 

+ -^r^(-A + r)g^+ -^r^(3 + 4r) ^ + — ^ 
6480 ^ ^ J ^^^^ g4gQ V ^ ; b/b3 ^ g^g 



1 1 

72 



1620 3240 
^ -r^(-5 + 24r + 24r2) —r^ Q'l 



3240 ' ' 3240 
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1 1 1 1 

80 



3240 486 2160 927 

5.5. Magnetic moments and susceptibility 

To obtain the magnetic moments /i„ we can use a special method which is available if 
one knows the moments t„ for all values of the coupling constants J^jjy. We replace H 
by the one parameter family of Hamiltonians = H + ^ (S^ — Nr). Equivalently we 
can substitute J^^ i— )■ J^j, + a for all coupling constants. The magnetic moments then 
result from differentiating Tr(if""'"^) w. r. t. a and finally setting a = 0: 



r) -I- 1 

= '-^Ti (HSiS' - Nr)) (46) 

°^^ <""^"f"^'" (3.,.-iV.U. (47) 

We can calculate the left hand side of ( H6|l if we insert the results for the moments 
and consider "derivatives" Q' of multigraphs defined in the following way. Let g^'^^ denote 
the multigraph Q but with one bond removed, J\f{i,j) h- )■ Af(i,j) — 1. If Af{i,j) = 
then we set Q^'^^^ = 0. Further let G{Q) and G{Q^''^^) denote the respective symmetry 
groups. Then we define 

One has, so to speak, to break each bond of the multigraph and to sum over all results. 
Further, one has to introduce factors which compensate for the possible change of sym- 
metries. For example, A ' = 6 A + It is obvious that the evaluation of Q' just 
yields ^ ^|q-o' Then it is a straightforward task to calculate the magnetic moments 
/io, • • • ,/^5 by using the above results for the tn- We will not, however, display these 
results and pass to the c„. 

The coefficients of the high temperature expansion of x = P '^'^^Tricy.^-^pH))^'^ '"^^ 
be expressed through the and the tn which occur as coefficients of the series in the 
numerator or in the denominator, respectively. The first 6 coefficients are given by: 

oo 

X=5^c„/3'^ (49) 

n=l 

= /io /3 - Ail + ^(/i2 - /io h)^'^ + ^(^sAto + 3^2/^1 - /is)/^'^ 
z 6 

+ ^(6^2Aio - U^^o - 4^3/^1 - 6t2Ai2 + /i4)/3^ 

+ Y^(^5/io - SOt^/ii + 5t4/il + lOts/is + 10t2(-2t3/io + /is) - /i5)/3^(50) 
+ .... 
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Inserting the known values for the t„ and the yU.„ yields the desired results for the 
c„. Similarly as in section 13. 2[ a variety of product rules can be used to simplify the 
resulting expressions revealing the extensive character of the c^. By this we mean the 
following. If the spin system under consideration would have a periodic lattice structure 
of, say, K lattice units with periodic boundary conditions, it follows immediately that 
the evaluation of a single multigraph Q linearly scales with and hence with A^, as 
long as Q is connected. For unconnected Q the evaluation scales with K'^ where c is the 
number of connected components of Q. Obviously, products of evaluations of connected 
multigraphs Qfj, would scale with K"^. It turns out that the elimination of these and 
higher products in the expression for the c„ by means of the rules (HI)- (!2T!l also eliminates 
the evaluation terms of unconnected multigraphs. This has to be expected on physical 
grounds, since the total susceptibility of a spin lattice should be an extensive quantity, 
i. e. linearly scale with K. But it is an additional consistency test of our results that 
the non-extensive contributions to the c„ actually cancel. 

We will now represent the results for the susceptibility's HTE up to sixth order in 
the inverse temperature /3. ci and C2 are well-known, C3 has already been published in 
[3], but C4,C5 and Cg seem to be new, although the polynomials which appear in these 
expressions are known up to 8th order, see [1]. 

Nr , , 

Cl = ^, (51) 

C2= (52) 

C4 = Y^^'(-l + r + r^) ^ + Ir^ ^ + ^r^ ^ - ^r"^ ^, (54) 

^5 = 7I7^^'(-3 + 6r + 4r2) ^ - ^r\-\ + r + r^) ^ 
o48 4U5 

+ —r'-^Gu —r^Gi3 ^r^(-3 + 8r + 8r2)^ 

108 243 ^ 540 ^ ' 

^ r^^--i-r^^--^r^^--i-r^^+-^r5^, (55) 



486 162 243 243 243 

23 



C6 = — r^(192 - 522r - 67r^ + 240r^ + 160r^) Q. 

113400 ^ ^ 



X 1 _ 

-r^(-3 + 6r + Ar"^) ^24 - 7^r-\-3 + 3r + 2r^) G: 



1944 ' ' 972 

G27 r^{-Q + 21r + 16r^) 



26 



972 2430 
1 9 



28 



+ ^— r^(-3 + 8r + 8r^) ^29 + — ^r*(-l + r + r') ^30 
4860 ^ ' 1215 ^ ' 

^ -r^G^+ —^r%9 - 126r - 12r^ + 152r^)^- — 



324' ' 24300' ' ' ' ■ " ■ ' — 648' 
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1 A / „ . „ „ On — 5 /I — 1 ^ ^ ^ 

39 



1 2 1 

486 729 1620 ^ ^ 
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